An all silicon quantum computer 
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A solid-state implementation of a quantum computer composed entirely of silicon is proposed. 
Qubits are Si nuclear spins arranged as chains in a 28 Si (spin-0) matrix with Larmor frequencies 
separated by a large magnetic field gradient. No impurity dopants or electrical contacts are needed. 
Initialization is accomplished by optical pumping, algorithmic cooling, and pseudo-pure state tech- 
niques. Magnetic resonance force microscopy is used for readout. This proposal takes advantage 
of many of the successful aspects of solution NMR quantum computation, including ensemble mea- 
surement, RF control, and long decoherence times, but it allows for more qubits and improved 
initialization. 
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The primary difficulty in the construction of quantum 
computers is the need to isolate the qubits from the 
environment to prevent decoherence, while still allow- 
ing initialization, control, and measurement. To date, 
the most successful experimental realizations of multi- 
qubit, many-gate quantum computers have used room- 
temperature, liquid NMR with "pseudo-pure" states [Q. 
These computers are able to maintain isolation from the 
control and measurement circuitry by employing weak 
measurement on a large ensemble of ~10 18 uncoupled, 
identical molecules. Although such a large, highly mixed 
ensemble may bring the existence of entanglement into 
question , the arbitrary unitary evolution afforded by 
the RF-controlled quantum gates assures that these com- 
puters behave non-classically Their principal limi- 
tation results from their small initial nuclear polariza- 
tion. The size of the effective sub-ensemble of nuclei 
contributing to the pseudo-pure state, and hence the ef- 
fective Signal-to-Noise Ratio (SNR), decreases exponen- 
tially with each added qubit, leaving this method unlikely 
to exceed the 10-qubit level without substantial modifi- 
cation j|. 

The proposals of Kane |5| and others to use single nu- 
clear spins in a low-temperature solid solve the scalability 
problem of solution NMR, but they introduce the prob- 
lem of single-nuclear-spin measurement. It remains an 
experimental challenge to fabricate a structure in which 
individual nuclei are sufficiently coupled to an electronic 
system for single-spin measurement, but also sufficiently 
isolated for long coherence times. 

In this Letter, we propose a different solid-state NMR 
implementation of quantum computation which intro- 
duces electron-mediated cooling, but maintains the weak 
ensemble measurement that has made solution NMR 
quantum computers so successful. The device is made 
entirely of silicon, with no electrical gates or impurities. 
As will be discussed below, the qubits are spin-1/2 nuclei 



that are located in relatively isolated atomic chains, as 
shown in Fig. [l[ The nuclei within each chain are distin- 
guished by a large magnetic field gradient created with 
a nearby microfabricated ferromagnet Each nucleus 
has about 10 5 ensemble copies in a plane orthogonal to its 
chain. This structure is embedded in a thin bridge whose 
oscillations provide readout via magnetic resonance force 
microscopy (MRFM) 0, §. 

The advantages of an all silicon implementation are 
many. Foremost, the crystal growth and processing tech- 
nology for silicon are highly matured. In particular, the 
most sensitive structures for force detection have been 
made from pure silicon. Also, it is good fortune that the 
family of stable nuclear isotopes of silicon is quite simple: 
95.33% of natural silicon is 28 Si or 30 Si, which are both 
spin-0, and 4.67% is 29 Si, which is spin-1/2, perfect for 
the qubit. Thus, silicon is well suited for nuclear-spin iso- 
tope engineering Another crucial motivation for the 
choice of silicon is the observation that nuclei in a semi- 
conductor may be polarized by cross-relaxation with op- 
tically excited, spin-polarized conduction electrons [fl0| . 
Although there are many other means for dynamic nu- 
clear polarization, optical pumping in semiconductors 
has one important feature: the electrons whose hyper- 
fine couplings make the polarization possible recombine 
shortly thereafter and, hence, do not contribute to de- 
coherence during the computation. The absence of any 
spurious spins, nuclear or electronic, should leave the 29 Si 
nuclei well decoupled from their environment. 

In the following, we describe an example procedure for 
fabricating the structures of Fig. [j] |y]] . We start with an 
isotopically depleted 28 Si(lll) (< 1% 29 Si) wafer which 
has been miscut 1° towards (112). Oxygen atoms are 
ion-implanted from the surface of the substrate followed 
by appropriate annealing in order to form a buried oxide 
layer, which will be removed later to form an open space 
below the vibrating bridge. The wafer is transferred to a 



2 



Dy micromagnet 




FIG. 1: The figure shows the integrated micromagnet and 
bridge structure. The bridge has length I = 300 /im, width 
w — 4 /im, and thickness t = 0.25 /im. The micromagnet 
has length L — 400/im, width W = 4 (jm, and height // = 
10 /im, and produces a field gradient of dB z jdz — 1.4 T/fim, 
uniform over a 100 /im by 0.2 /im region inside the bridge. 
The insert shows the structure of the silicon matrix and the 
terrace edge. The darkened spheres represent the 29 Si nuclei, 
which preferentially bind at the edge of the Si step. 



multi-chamber molecular beam epitaxy (MBE) machine 
equipped with a scanning tunneling microscope (STM) 
and a pre-heating stage. It has been demonstrated that 
highly regular arrays of steps can be produced on vici- 
nal Si(lll)7x7 using a simple multi-step annealing se- 
quence |12|. Due to the high energy of interrupting or 
misplacing 7x7 domains, this procedure leads to atomi- 
cally straight step-edges along the (110) direction for up 
to 2 x 10 4 lattice sites. The average terrace width is 
~15 nm. Once we verify the straight step edges using 
the STM, we transfer the wafer to the growth chamber 
to form atomic chains of 29 Si using the step-flow mode at 
temperatures T ~ 850°C. Arrival of 29 Si isotopes at the 
substrate surface induces a surface transition from 7x7 



to lxl i.e., 29 Si isotopes travel to the step edges 
and form atomically straight lines. We terminate the 
evaporation of 29 Si when the atomic chains are one atom 
wide and run completely along the step. A 28 Si capping 
layer of 15 nm is grown on top of 29 Si chains and we re- 
peat the same multi-step annealing, 29 Si chain growth, 
and capping sequence to produce replicas of parallel 29 Si 
chains. The last step in the growth process is the capping 
of ensembles of 29 Si chains by a thick 28 Si layer. 

Once formation of the 29 Si wire block is complete, the 
narrow bridge is created with e-beam lithography. A 



plasma etcher removes unprotected silicon down to the 
buried oxide layer. An HF vapor etch or acid solution 
removes the exposed oxide to release the structure; it is 
then dried in a critical point dryer. After etching 3 /im 
further into the substrate and through the oxide layer, 
the dysprosium (Dy) micromagnet can be evaporatively 
deposited s — 2.1 /im from the bridge and defined as a 
parallelopiped using a shadow mask. 

Calculations similar to those in Ref. || show that 
the magnetic field gradient due to the micromagnet is 
dB z /dz — 1.4 T//im, which is persistent over the thick- 
ness of the bridge and is superposed with a large homo- 
geneous field Bq of ~ 7 T. The distance in the z-direction 
between two 29 Si nuclei in an atomic chain, which we no- 
tate a, is 1.9 A, so the gradient leads to a qubit-qubit fre- 
quency difference of Aus — a-fdB z /dz = 2tt x 2 kHz. The 
active region is a 100 /im by 0.2 /xm area in the center of 
the bridge, containing N = 10 5 chains persisting over the 
bridge thickness. This active region is sufficiently small 
and the magnetic field sufficiently homogeneous that all 
equivalent qubits in an atomic plane lie within a band- 
width of 0.6 kHz. 

For initialization of the quantum computer, we pro- 
pose to employ optical pumping, algorithmic cooling, and 
pseudo-pure state techniques. The premise of optical 
pumping is that nuclei exchange Zeeman energy with a 
bath of electrons which have been preferentially excited 
into a single spin state by circularly polarized light. The 
nuclei thereby relax thermally to an effective spin tem- 
perature corresponding to the non-equilibrium electron- 
spin polarization. Once those electrons recombine, the 
nuclei retain their spin polarization for the "dark" T\ 
time, which is extremely long (200 hours in Ref. fL0[). 
In low-field (~ 1 G) experiments at 77 K Jh], [jdj , nuclear 
polarizations have not exceeded 0.1% due to limited elec- 
tron spin polarizations and long recombination times, a 
result of silicon's indirect bandgap. Improved nuclear 
polarization in silicon may be observable in higher mag- 
netic fields (~10 T) and lower temperatures (~1 K), and 
in silicon nanostructures where rapid recombination of 
electrons via surface states can help maintain the elec- 
tron spin polarization Q . 

The physical cooling afforded by optical pumping is 
well suited to complement an algorithmic cooling tech- 
nique introduced by Schulman and Vazirani jl^, [lTj . This 
technique redistributes the entropy among a register of 
qubits to a subregister that is then discarded (decoupled 
and ignored) [fijl. This method is expensive; it takes time 
to perform the (classical) logic operations and, worse, it 
sacrifices many qubits. An initial register of size no will, 
for small initial polarization po, shrink to size uqPq/ 2 hi2 
if the procedure is taken to the entropy limit. However, 
the very long T\ affords ample time for the procedure, 
and the number of available qubits in our configuration 
before initialization can be thousands. Moreover, the al- 
gorithm need not be taken to the entropy limit, since 
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large but still incomplete polarizations can be handled 
with pseudo-pure state techniques, the consequences of 
which will be discussed below in the context of SNR. 

The secular component of the dipolar Hamiltonian 
which couples the ith spin to the jth spin within one 
chain is written 11191 
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where r.y is the length of the vector connecting the spins 
and 9ij is its angle with the applied field. Nearest neigh- 
bors along the proposed atomic chains are not exactly 
parallel to (110), but rather zig-zag with angle 
satisfying cos 2 0i,i+i = 2/3, leaving 5uj = = 
2ir x 0.4 kHz. Other terms of the dipolar Hamiltonian 
require the exchange of energy in the amount of HAlu or 
more. The long T\ in silicon indicates the inefficiency by 
which this energy may be exchanged with degrees of free- 
dom external to the nuclei. In strongly coupled dipolar 
systems, this energy can be compensated for by the dipo- 
lar bath of identical nuclei |Q . In the present scheme, the 
members of the nuclear ensemble are so far removed that 
this energy bath is absent. Hence, we assume these non- 
secular terms are well suppressed. The Hamiltonian of 
Eq. (|l|) may be "switched off' by applying a periodic 
succession of narrow band it pulses at, for instance, the 
«th resonant frequency plfl . Simultaneous decoupling of 
more than two qubits may be accomplished by timing 
the selective n pulses according to the entries of an ap- 
propriately sized Hadamard matrix; a pair of qubits may 
be selectively recoupled in order to implement two-bit 
gates . Note that this Hadamard pulse scheme serves 
the second purpose of refocussing inhomogeneous broad- 
ening caused by the in-plane nonuniformity of the field 
gradient and any small bulk susceptibility effects. 

The presence of a large magnetic field gradient pro- 
vides a natural means for performing MRFM on a mag- 
netization M z , since this technique is sensitive to the 
gradient force given by F z — M z dB z / dz. The experi- 
ment is performed in high vacuum (< 10 -5 torr) and at 
low temperatures (4 K). A coil is used to generate the 
RF pulses for logic operations and decoupling sequences; 
it also generates the continuous-wave radiation for read- 
out. An optical fiber-based displacement sensor is used 
to monitor deflection of the bridge using interferometry. 
Sub- Angstrom oscillations can be detected; larger oscil- 
lations can be damped with active feedback which avoids 
additional broadening while maintaining high scnsitiv- 
ity[f§. 

Readout is performed using cyclic adiabatic inversion 
p9[ , which modulates the magnetization of a plane of nu- 
clei at a frequency near or on resonance with the bridge. 
The spins of resonant frequency are irradiated with the 
RF field B x = 2BiCOs{uJit — (f2/u; m ) cos(u m t)}, where 



uj m is the modulation frequency chosen to be near the 
resonance of bridge oscillations, and f2 is the frequency 
excursion, which should be much smaller than Aoj ||. 
The z component of the ith plane's magnetization is de- 
duced from the phase of the resulting bridge oscillation. 
Simultaneous detection of signals from multiple planes is 
possible if the different planes to be measured are driven 
at distinct modulation frequencies u> m . 

The force resolution for MRFM is limited by thermal 
fluctuations of the mechanical oscillator pifl . Force res- 
olutions of 5.6 x 10 -18 N/VHz have been reported for 
single crystal silicon cantilevers at 4 K j2f| . The thermal 
noise threshold of the bridge structure in Fig. |] is esti- 
mated to be ~ 1.2 x 10~ 17 N/VHz based on a lumped 
harmonic oscillator model and assuming a modest qual- 
ity factor Q of 10 4 j26). This model yields an estimated 
spring constant of k i=s 0.0042 N/m and a resonance fre- 
quency lu c /2tt 23 kHz. The detectable signal will de- 
pend upon the initial polarization p after optical pumping 
and algorithmic cooling. The force from the subensem- 
ble magnetization corresponding to a pseudo-pure state 
is estimated (3, M as 
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The number of qubits available in this scheme may be 
found by maximizing n in Eq. (Q) such that the force 
exceeds the thermal noise threshold; the results of such 
maximization are shown in Fig. ^. At low p, exponential 
improvements in p are needed to increase the number of 
measurable qubits n. Once p exceeds about 60%, how- 
ever, we find n ~ (l+p)/(l — p), escaping the exponential 
downscaling which plagues solution NMR. 

There are several possible sources of decoherence in 
this proposal: magnetic fluctuations in the dysprosium, 
thermal currents in the dysprosium, fluctuations of para- 
magnetic impurities in the silicon, thermal motion of the 
bridge, and uncontrolled dipolar couplings between nu- 
clei. We estimate that the latter two sources are the most 
important, so we limit our present discussion to them. 

A calculation using the thermal noise statistics of a 
high-sensitivity mechanical oscillator |24|] estimates the 
Ti timescale due to the bridge's thermal drift as T| = 
kuj c Qa 2 / Auj 2 k-QT ss 25 s. Active feedback stabilization 
is expected to increase this timescale by as many as four 
orders of magnitude, bringing it close to the lengthy T\ 
timescale. 

There are two sorts of uncontrolled dipolar couplings 
to consider. There are the in-plane couplings, in which 
the participating nuclei have equal Larmor frequencies 
and can thus participate in spin flip-flop processes. Since 
the average distance between chains is ^15 nm, these 
processes cause decoherence on a time scale of order 
T 2 h - 4tt(15 nm) 3 /7 2 fyio - 100 s. This already lengthy 
timescale may be increased by the addition of dipolar 
refocussing sequences such as WAHUHA tH] . 
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FIG. 2: A plot showing the scalability of the present scheme. 
The solid curve, corresponding to the left axis, shows the 
polarization p needed in order for a number of qubits n to be 
measurable. The dashed lines, corresponding to the right axis, 
plot the number of logic gates (decoherence time T2 divided 
by the pulse sequence cycle time t c ) times the length of a 
"block" of the decoupling sequence L against n for several 
values of T§. The inset, showing F(m), is explained in the 
text. 



There will also be spurious couplings between nuclei in 
one homogeneous plane and copies of nuclei in the next. 
Such couplings occur when two qubits are recoupled for a 
logic gate, and they cause a small error in each logic gate. 
To estimate this error, suppose we recouple one qubit to 
another qubit to planes away. The strongest coupling 
seen by a nucleus is the in-chain coupling, and its rate is 
~ 5uj/m 3 . The couplings to all of the neighbors may be 
treated as a T2 decoherence process, with 
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Here, T2 has been estimated as the inverse square root 
of the second moment [^Tj and A,; is the lateral distance 
to the ith chain normalized by a. The ratio of the ap- 
proximate gate time m 3 /Soj to T\ m is plotted in the in- 
set of Fig. H This function, F(m), represents the ap- 
proximate error in the gate when attempting to couple 
qubits to planes apart. To keep gate errors low, the com- 
puter should couple only nearby neighbors and handle 
more distant couplings by bit swapping. In this way, 
gate errors due to unrefocussed nuclear couplings be- 
tween chains is limited to approximately F(l) « 10 -6 . 

These cross-couplings also influence the clock speed of 
this computer, as we now explain. When qubits are de- 
coupled by the Hadamard scheme mentioned above, the 
resulting pulse sequence has a clock time t c = Ln 2 /Auj, 
where n is the number of qubits being decoupled and 
L/Aw is the amount of time devoted to one n pulse [ p2| . 
As qubits are added, the amount of time needed to recou- 
ple them grows. For very large n, however, some qubits in 
the chain become so distant that the rfj factor in Eq. (|l|) 
renders their interaction negligibly small. In this case 



the Hadamard pulse scheme can be truncated to decou- 
ple qubits only in sets of /, and a new decoherence process 
is introduced with T^Sco = l 3 [l + F 2 (7)]~ 1/2 . If we com- 
pare the total decoherence time T 2 = (1/T§ + l/T^) -1 , 
where T° combines all other decoherence processes, to 
the clock speed t c , we find that there is an optimum I 
at which to truncate the Hadamard pulse scheme. The 
effective number of logic gates T 2 /t c , therefore, at first 
decreases as n 2 and then flattens once n reaches this op- 
timum I, as shown in Fig. |[ 

By consulting Fig. [| we see that this scheme can, for 
sufficiently large T® and sufficiently high polarizations, 
allow substantially more qubits than solution NMR, 
without the need for single-spin measurement or unre- 
alistic advances in fabrication, measurement, or control 
technologies. 
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